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Abstract 

Integrable defects in two-dimensional integrable models are purely transmitting thus 
topological. By fusing them to integrable boundaries new integrable boundary conditions 
can be generated, and, from the comparison of the two solved boundary theories, explicit 
solutions of defect models can be extracted. This idea is used to determine the trans- 
mission factors and defect energies of topological defects in sinh-Gordon and Lee- Yang 
models. The transmission factors are checked in Lagrangian perturbation theory in the 
sinh-Gordon case, while the defect energies are checked against defect thermodynamic 
Bethe ansatz equations derived to describe the ground-state energy of diagonal defect 
systems on a cylinder. Defect bootstrap equations are also analyzed and are closed by 
determining the spectrum of defect bound-states in the Lee- Yang model. 

1 Introduction 

Recently, there has been an increasing interest in integrable quantum field theories including 
defects or impurities. This is motivated both by the realistic physical applications in statis- 
tical and solid state physics and also by the need of theoretical understanding of this so-far 
unexplored field. 

The community of integrable systems have not payed much attention to defect theories at 
the beginning due to the no-go theorem formulated by Delfino, Mussardo and Simonetti in 
[HE]. The theorem, formulated originally for diagonal theories and extended later for a large 
class of non-diagonal ones in [3], states that a relativistically invariant theory with a non-free 
integrable interaction in the bulk can allow only two types of integrable defects: the purely 
reflecting and the purely transmitting ones. (Although some effort has been made to overcome 
this obstacle by giving up Lorentz invariance, see for instance [I] and references therein, in 
the present paper we restrict ourselves to the relativistically invariant case.) 

The analysis of boundary integrable theories was initiated in [5] by formulating, in an 
axiomatic way, the properties of the reflection matrix: unitarity, boundary crossing unitarity 
and boundary bootstrap equation. The boundary bootstrap framework was completed by 
introducing boundary Coleman-Thun mechanism [6] and the bulk bootstrap equations [7|. 
Later this framework got a sound basis by developing boundary quantum field theories from 
first principles in [HI [9]. The success of the boundary bootstrap approach resulted in a large 



class of closed bootstrap theories in which the boundary reflection factors together with the 
spectrum of boundary excited states were determined [5J El UHl ttH El EH E] • The solutions, 
obtained by the bootstrap method, are not connected, however, to other formulations of 
the model such as to the classical field theory or to the perturbed conformal field theory 
(both defined by a Lagrangian). To connect the different descriptions one either has to check 
the reflection factors perturbatively, like in [15], or solve the theories in finite volume. The 
boundary thermodynamic Bethe ansatz (BTBA), developed in [16], systematically sums up 
the finite size corrections by taking into account the scatterings and reflections. By analyzing 
its small volume limit the needed link between the bootstrap and perturbed conformal field 
theoretical descriptions can be established. 

As the no-go theorem showed non-free integrable defect theories are purely transmitting. 
This fact kept back the researchers for some time to analyze these models until new life was 
put into the subject due to their explicit Lagrangian realizations [T7]. Following the original 
idea many integrable defect theories were constructed at the classical level [IB], [19|, [20] . The 
basis for the quantum formulation of defect theories is provided by the folding trick [21] by 
which one can map any defect theory into a boundary one. As a consequence defect unitarity, 
defect crossing symmetry and defect bootstrap equations together with defect Coleman-Thun 
mechanism are derived. Despite of these results the explicitly solved relativistically invariant 
defect quantum field theories are quite rare, containing basically the sine-Gordon and affine 
Toda field theories [22], [23] El] and even in these cases the explicit relation to their Lagrangian 
have not been worked out yet. 

One may think that purely transmitting theories are too simple and there is no point 
to analyze them, but we would like to argue that they carry very important information 
about an integrable quantum field theory, without which our knowledge cannot be complete. 
Purely transitivity implies the conservation of momentum from which the topological nature 
of the defect follows. Thus, such defects can freely be transported in space without affecting 
the physics of the theory. We can either move them close to each other or move them to 
integrable boundary conditions and, as a result, new integrable boundary conditions can be 
generated. These ideas were successfully applied in conformal field theories [25] [26] [2?], in 
integrable lattice models [28] [29] and the aim of the present paper is to exploit it in solving 
integrable defects in the sinh-Gordon and Lee- Yang theories. 

The paper is organized as follows: In section 2, on the example of the sinh-Gordon the- 
ory, we show how new boundary conditions can be obtained by fusing integrable defects to 
boundaries at the classical level. We focus on two cases in detail: fusing the integrable defect 
to Dirichlet boundary condition (DBC) the perturbed Neumann boundary condition (PNBC) 
can be obtained, while fusing it to the Neumann BC a new class of time-dependent integrable 
BCs can be generated. Section 3 recalls the quantum version of the fusion method together 
with the properties of transmission factors. By comparing the already known reflection factors 
and boundary energies of the DBC to those of the PNBC solutions for the defect transmission 
factor and defect energy can be extracted. The same method is then used to determine defect 
energies and transmission factors of the scaling Lee- Yang model. This method not only pro- 
vides the explicit solutions of the sinh-Gordon defect theory but also relates its parameter to 
that of the Lagrangian. Since the relation obtained here is different from the suggestion of [23] 
we perform a perturbative analysis at one-loop level in section 4. (In the subsequent paper [21"] 
the same authors raised the possibility of the quantum renormalization of the transmission pa- 
rameter, which is confirmed at one-loop level here). In section 5 the calculated defect energies 
are subject to another consistency check. For this we derive a TBA equation to describe the 
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groundstate energy of a diagonal defect system on a cylinder. Prom a careful ultra-violet (UV) 
analysis defect energies are extracted and the previous results are verified. By completing the 
defect bootstrap program we analyze the singularity structure of the transmission factors in 
section 6. Since the sinh-Gordon transmission factor does not contain any singularity in the 
physical strip we analyze the Lee- Yang model only. For each pole of the transmission factor in 
the physical strip we associate either a defect boundstate or a defect Coleman-Thun diagram 
and calculate the excited transmission factors in the former case from the defect bootstrap 
equation. Finally, we conclude in section 7 and give directions for further research. 

2 Fusion method at the Lagrangian level 

In this section we demonstrate, on the example of the sinh-Gordon (ShG) theory, how new 
integrable boundary conditions can be obtained by the fusion method at the Lagrangian or 
classical level [30J. 

The ShG theory is defined on the whole line by the following Lagrangian 

1 1 Tfl 

W*) = - 2 (d&? - -{d x <S>? - cosh 6<D , (1) 

It describes an integrable field theory and by restricting the theory to the half line this property 
can be only maintained, whenever the following boundary potential is introduced [5] 

£ BS hG = &(-x)£ ShG (<f>) - 6(x)B($) ; B($) = M cosh~($-^ ) (2) 

By varying Mq from to oo the arising boundary condition interpolates between the Neumann 
d x Q\ x= o = and the Dirichlet $(0,i) = (fo ones. 

The most general integrable defect condition can be obtained by the analytical continuation 
of the sine-Gordon result [23J. The Lagrangian in the ShG case reads as 

£ DS hG = e(-x)£shG(*-) - *(x)Z>(*_,*+) + e(x)£ S hG($+) (3) 
where the defect potential contains just one single parameter: 

2 J D($_,$+) = - + M d e^ cosh ^($ + + +M d e"' i cosh ^(<£» + - 

Here & T are the fields living on the left/right half-line, respectively and M c \ = ^f 1 - 

The fusion idea is based on the integrability of the defect: Integrability guaranties the 

existence of an infinite number of commuting conserved charges which results in the possibility 

of shifting the trajectories of particles, without changing the amplitude of any scattering 

process. The shifting of all the trajectories can alternatively be described by shifting the 

location of the defect, which then, does not alter the physics. 

At the level of the Lagrangian this observation can be formulated in the following way: The 

spectrum of the system, which contains a defect in the origin, x = 0, in front of a boundary, 

located at x = a, 

£ D BShG = 8(-ar)£ ShG (*-) - 6(z)D($-,$ + ) + 0(ar)0(a - x)£ ShG ($ + ) - 8(x - a)B(<S> + ) 
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does not actually depend on a. Thus we can perform the a — > limit and represent the same 
system as a boundary one, but with a different (dressed) boundary condition: 

CnBShG = 0{-x)£shG($-) ~ S(x)B' (*_,*+) ; B' (*_,*+) =£>($_,$+)+£($+) 

The field $ + lives on the boundary only and can be thought naively to be a boundary degree 
of freedom. It does not have, however, any kinetic term so it merely implements a new 
time-dependent integrable boundary condition. Let us specify these findings in two concrete 
examples that will be used later on. 

If the original boundary condition is the Dirichlet one with $(a,i) = <po, then the arising 
dressed boundary condition is 

<&+) = £>(*- , <p ) = C ° sh + 0o) + °° sh ^ ($ - " 0o) (4) 

where we dropped the total time derivatives. This BC is exactly of the form of ([2j) with 
parameters 

M c icosh( M ±^o) = M e^ 

Thus by fusing the integrable defect to the DBC we can reconstruct the most general (two pa- 
rameter family of) PNBCs. Interestingly (fio and fj, are the classical analogues of the parameters 
in which the boundary reflection is factorized [45J, see also (|13)14p in section 3. 
By fusing the defect to the NBC we obtain the boundary potential 

B' (*_,*+) = £>($_,$+) 

Variation of action provides BCs in the form: 

dB'($+ $ ) dB'(®+ <E> ) 

**-\**> = - ; m + -d x *-)\ x=0 = a "^*- ) (5) 

By expressing <I> + in terms of $_ and dt&- then plugging back to the second equation we 
obtain a highly non-trivial boundary condition for containing its second time derivative, 
which is nevertheless integrable as it follows from the construction. Obviously, this solution 
was not covered by the two parameter family of (time-independent) integrable boundary con- 
ditions determined in [5], thus by the fusion method we were able to construct a new type 
of integrable BC. By fusing other integrable defects with new free parameters to this dressed 
boundary we can generate integrable BCs with as many parameters as we want. What is nice 
in the construction, that the solution of the defects transmission factor will provide, via the 
fusion method, solutions for these general integrable BCs, too, as we will show in the next 
section. 

Finally, we note that similar construction can be used in the case of the Lee- Yang model, 
however, the explicit form of the integrable defect perturbation has not been identified at the 
Lagrangian level yet ( for details see the next section). 

3 Fusion method in the bootstrap 

For simplicity we present the fusion idea in the case of an integrable diagonal scattering theory 
with one particle type of mass m. The general discussion can be found in |21j. 
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In integrable bulk theories multi-particle scattering processes factorize into the product of 
two particle scatterings, S(9i2), where 9±2 = 61 — 62 is the rapidity difference of the scattering 
particles whose momenta are parametrized as pi = msinh#j. In relativistically invariant 
theories the two particle S'-matrix satisfies unitarity and crossing symmetry 

S(-6) = S- 1 (6) ; S(i7r -6) = S(6) 

Once boundaries are introduced the basic process is the multi-particle reflection. Integra- 
bility ensures its factorization into pairwise scatterings S(9ij) and individual reflections R(9i), 
where 9i is the rapidity of the reflected particle. The reflection matrix satisfies unitarity and 
boundary crossing unitarity [5] 

R{ - 9 ) = R-\6) ■ rC^-6) = S(26)rC^ + 6) 

Integrable non-free defects are severely restricted: they are either purely reflecting (thus 
boundaries, like above) or purely transmitting. The latter case can be described by the two, 
left (— ) and right (+), transmission matrices T_{6) and T + (—6). We parametrize T + such a 
way that for its physical domain (6 < 0) its argument is always positive. Transmission factors 
satisfy unitarity and defect crossing symmetry [21] 

T + (-6) = TZ l (6) ; IL(0) = T+(mt - 6) (6) 

If we place a defect with transmission matrices T±(0) in front of a boundary with reflection 
matrix R{6) then the fused boundary will also be integrable and have reflection factor R (6): 

R'(6) =T + (6)R(6)T_(6) (7) 

The correspondence j7]) between the original R{6) and the fused R' (6) reflection factors can 
be used either to generate new BCs or, if the two BCs are already known, to solve defect 
transmission factors. This will be illustrated in the next subsections for the sinh-Gordon and 
Lee- Yang models. 



3.1 Solution of defect sinh-Gordon theory 

The spectrum of ShG theory defined by ([T|) consists of one particle type with mass [32] 



40F / -mlrjl + i 2 ) ^ b 2 

m - r (i^)r(i + |) V b 2 r(-b 2 ) J ; ^"8^ + 52 

The two particle scattering matrix is given by 



ITTX \ 



Binhe + ismBir sinh(| — ^) 

It has no poles in the physical strip: < 6 < in and is invariant under the weak-strong duality, 
— > |f . The bulk energy density turns out to be [33] 



Integrable boundary conditions can be either Dirichlet type with $(0, t) = 4>q or PN type 
(EJ. The corresponding reflection factors can be obtained from the analytical continuation of 
the sine-Gordon's first breather's one [5j [331 ■ ^ n l ne Dirichlet case it reads as 

where the reflection parameter r/Dir is related to cj>Q as 
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b 

The boundary energy has been also calculated [T6 



VDir = -r4>0 (11) 



Dir m ( nB TlB \ 

e bdry (VDir) = 4^5^ I 2 COsh Br] DiT — Sill — COS — 11 (12) 



In the PN case the reflection factor turns out to be 



lw l 



= ^TTTBf 7^ ;( hag ■ n (13) 

^ 2 2 ' ^ 7T + 2 7 

while the relation of ij, ?9 to the parameters of the Lagrangian (UV-IR relation) is 



b z b V0 / 2 

Mcosh— (t/±0) = M e T 2 ; M = m d J . (14) 
87r y o" 1 sm(tr/8) 

The boundary energy has been also determined [31, 45j as 



e™rv(^ i?) = , m „ (2 cosh #77 + 2 cosh Bi9 - sin ^— - cos ^— - 1 ) (15) 
' 4 sm i37r V 2 2 



We note that the results - both for the UV-IR relation and for the boundary energy - were 
obtained in the framework of perturbed BCFT in which the perturbing operator is normal- 
ordered to have a definite scaling dimension. 

The integrable defect potential for the sinh-Gordon model can be written as in Let us 
denote the transmission factors by T±(6,fj,). Fusing classically this defect to a DBC a PNBC 
can be obtained . The quantum analogue of this statement in view of (JJJ) is 

R PN (6, 7?, 0) = T+(0, fx)Rnn{e, r/ D ir)T_(#, fx) 

Comparing the reflection factor of the PNBC ffT3|) to that of the Dirichlet one (fTOj) and taking 
into account defect unitarity and defect crossing symmetry ([6]) we can extract the transmission 
factors for the defect. The simplest possible solution corresponds to rj = r/Dir and 



T-{6) = -i )| £ Ijf ; T+{0)=i )| * 1^ 16 

sinh(f + f + ^) +W sin h(f + f-^) 

All other solutions contain additional CDD type factors satisfying Actually the solution 
(fl~6j) itself is a CDD factor, therefore it is the simplest non-trivial solution of (EI). 
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To find the correspondence between the parameter of the quantum transmission factor i? 
and the Lagrangian parameter /x we follow the following strategy: Since the boundary results 
are derived in the perturbed BCFT normalization we allow not only the parameter \x but also 
M c \ to renormalize. Their renormalized quantum values are determined from the requirement 
that when fusing the defect to the DBC with (HU we obtain the PNBC with ([Til) . The unique 
solution turns out to be 

Me ±b i£ = M cl e^ (17) 

The renormalization of M c \ may depend on the scheme in which the quantum potential is 
defined. The b — > limit, in which M — > M c \, shows that $ is the quantum renormalized 
version of /J,. 

Also the defect energy can be extracted as the difference of the boundary energies corre- 
sponding to the PN (fT5|) and to the Dirichlet (fT2"T) one: 



/ q\ PN I o\ Dir / x mCOshSi? 

eDefW = e b Sy(7/,i?) - e bdry (7?) = 2gini?7r 



Summarizing, by the fusion method we were able to solve the defect theory defined by the 
Lagrangian (J3j) : The transmission factors are (fT6l) . the defect energy is (fT8j) . and the bootstrap 
parameter ■& parametrizes the Lagrangian as (fTTl) . We spend the next section to provide 
consistency checks of this solution. 

Once the defect theory is solved we can use it to generate new integrable BCs from known 
ones. In the example presented in section[2]the defect with parameter fi was fused to the NBC 
to generate a more general integrable BC (EJ. The quantum version of this fusion dresses up 
the Neumann reflection factor 

n (a\ _ (2) (j; ~ 2") (2 ~ 2") / 10 \ 
RmW) - (|_|) (|"+f) 

to the reflection factor 

(1) (1 _ B\ n _ B\ (im _ i\ 
R(9,$)=T + (e,^))Rv(9)T„(9,^)) = ^ B ' J jj ^ J (20) 

(i-f) (l + f)(^ + 5) 

Thus we solved the more general integrable BC defined by J5]) without doing any serious 
calculation. It is important to note, that the extra factor appearing in (|20l) compared to (fT9j) 
is a CDD factor. Consequently, we have determined the physical meaning of the CDD factors 
appearing in the reflection factors: they represent integrable defects of the form of (J5]) standing 
in front of integrable boundaries. In principle, we can fuse as many integrable defects with 
various parameters as we want, the resulting theory can be solved and its reflection factor 
contains the corresponding boundary CDD factors. 

Finally, we note that placing two defects with parameters $± = ±i(l — ^) after each other 
both the scattering matrix and the energy of a standing particle can be reproduced. Thus the 
defect with imaginary parameter can be considered as a 'half particle. Similar phenomena 
was observed in [23J at the classical level. 

3.2 Solution of defect scaling Lee- Yang model 

The scaling Lee- Yang model can be defined as the perturbation of the -M(2,5) conformal min- 
imal model with central charge c = — It contains two modules of the Virasoro algebra 
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corresponding to the Id and the tp(z,z) primary fields with weight (0,0) and (— -|), re- 
spectively. The only relevant perturbation by the field (p results in the simplest scattering 
theory with one neutral particle of mass m and scattering matrix [35] 

sinh 9 + i sin £ fl\ {2 

S{9) - 



sinh — i sin ^ \ 3 J \ 3 
The pole at 9 = ^ shows that the particle can form a bound-state. The relation 



S{e + i\)S{9-£) = S{9) 



however, implies that the bound-state is the original particle itself and the bulk bootstrap is 
closed. The bulk energy constant is given by ebuik = ~J~^ m2 - 

We can impose two conformal invariant boundary conditions in the model (36J E7] . They 
can be labeled by I and $ and correspond to the highest weight representations of a single 
copy of the Virasoro algebra with weight and — g, respectively. Introducing the integrable 
bulk perturbations with the I conformal invariant boundary condition the integrability is 
maintained and the reflection factor of the particle can be written as 

while the boundary energy is given by e bdry = T ~ ■"■) • ^ * ne conformal invariant bound- 
ary condition corresponds to $ then additionally to the bulk perturbation we can introduce a 
one-parameter family of integrable boundary perturbations and the corresponding reflection 
factor turns out to be 



2) \6J \ 3J\6J\6J\6J\ 6 

while the boundary energy is ebdry = y (v^ — 1 + 2 sin ^). The boundary bound-states were 
analyzed in [6] where the boundary bootstrap program was carried out. 

The Lee- Yang model has two types of conformal defects [38], but only one of them admits 
relevant chiral defect fields. They have weights (— |, 0), (0, — g). We conjecture that perturbing 
in the bulk and with a certain combination of these defect fields we can maintain integrability 
and arrive at a purely transmitting theory. We plan to analyze this issue systematically in 
a forthcoming publication. Let us denote the transmission factors of this integrable defect 
by T±(9). Using that the fusion of the defect to the I boundary results in the perturbed $ 
boundary we have 

R b (9)=T + {9,b)R 1 {9)T-(9,b) 

This is supported by the fact that fusing the conformal defect to the I boundary we obtain 
the $ boundary. Since the particle appears as a bound-state in the two particle scattering 
process the transmission matrix satisfies the defect bootstrap equation [3]: 

T_(9+ l ^)T^9-^) = T^9) (21) 

Using this relation together with the defect unitarity and defect crossing symmetry © we can 
fix the transmission factor as 

T(9) = [b+l][b-l] ; [x]=i sinh d + i ff) (22 ) 

v ' 1 JL J L J sinh(| + iff -q) y ' 
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(Actually the inverse of the solution is also a solution but the two are related by the 6^6 + 6 
transformation). The defect energy, as in the sinh-Gordon case, can be obtained as 

I ■ to 

£def — Cbdry ~ ^bdry — msm — 

We are going to recover this expression from the UV analysis of defect TBA in section 5. We 
also note that the defect with parameter b = 3 behaves as a standing particle both from the 
energy and from the scattering point of view. 



4 Peturbative calculations 

In this section we check the exact solution of the ShG defect system in the free/classical 
(b — ► 0) limit, and develop a systematic perturbative expansion. The parameter b 2 plays the 
same role as H which can be seen by scaling it out from the Lagrangian via <3? — > M>. Since 
the classical groundstate <3? = is invariant under this scaling, the b 2 — ► limit corresponds 
both to the free and also to the classical limit. Moreover, the perturbative expansion in b 2 is 
equivalent both to the loop expansion and to the semi-classical approximation. 



4.1 Classical/free limit 

As a first step we identify the 6^0 limit of the defect Lagrangian j3|) as: 



C 



9(-ar) 
S(x) 



~(W 2 -^* 2 _ 



+ e(x) 



+ m c \ [cosh fji + $i) +2 sinh /i 



Then we expand the fields on the two sides of the defect in terms of creation/annihilation 
operators 



*±0M) 



dk 1 



a±(k)e 



ikx—iuj(k)t _|_ + 



+ a%{k)e 



-ikx+iu)(k)t 



u{k) 



2vr 2u(k) 

where the a, a + operators are adjoint of each other with commutators: 

[a±(k), aj(fe')] = 2-K2u(k)5(k - k') 
Imposing the defect condition (obtained by varying the action) at the origin 

±d t $± T d x ® T = rod (sinh \x $± + cosh// 
we can connect the creation/annihilation operators as 

m c i sinh /j, =F ioo(k) 



k 2 + 



m- 



a±(±k) = T T (k)a T (±k) ; T T (k) 



m c \ cosh fx — ik 



k > 



This shows that the defect is purely transmitting, that is we do not have any reflected wave. 
The transmission factor in the rapidity parametrization (k = m c \ sinh 9) can be written also 
in the following form: 

,sinh(f-f + f) 



T-(9) 



sinh(f + f + £ 
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Clearly it has exactly the same form as the exact quantum one (fT6l) except the Bfl <-> \x 
replacement. Having observed this coincidence the authors in [23] suggested that they might 
be the same B-d = fx. Using our defect UV-IR relation fjTT|) we can perform the expansion: 



« = rti- s + .. 



(23) 



The term of first order shows that our exact solution is correct in the classical limit, i.e. for 
b — > 0. The term of second order shows that the B$ = \x relation suggested in [2B] is not valid: 
B'd acquires nontrivial quantum correction. Since the renormalization of the parameter fi is 
crucial to decide about the two proposals we perform a perturbative check at order b 2 . 

4.2 Perturbation theory 

As a first step we collect the free propagators. If the fields are on the same side of the defect 
we have 



o(0\T^ ± (x,t)^(x',t)) |0) = J^-j-J^e**™') = C%M) 



where q = {k,u) and y = (x,t). The absence of an e lk ^ x+x ) term shows the absence of 
reflection. The other two point functions are 



o(0\T^(x,t)^ ± (x',t)) \0) = J^-^^T ± (u;,k)e^y-y^ = G^ y ,y') 



where 



T±(u,k) 



m c \ sinh fxrtiu 



m c \ cosh ii — ik 

In the final equations we used the fact that the u> contour can be closed on the upper/lower half 
plane. As it was shown in [8] the reflection factor can be read off from the {T($±$±)} propagator 
of the fields. The defect/boundary equivalence [21] then implies that the transmission factor 
can be read off from the (T(& T &±)} propagator. 
The perturbation at order b 2 follows from (pQ): 



SC. 



-e(-c) 



2 12 



mi-ib 



4! 



@(C) 



4! + 



-8(0 



m c \b 2 
4-4! 



[cosh fj, ($4 + 6®1<S> 2 - + $1) +4sinh^$+$_ + 



We calculate the propagators upto first order in b 2 : 

(0\T($ T (x,t)$±(x',t'j) |0) = (Q\T[*T:(x,t)*±(x,t)(l-i I d£ I drSC + ...)\ |0) 



Using Wick's theorem we obtain the contribution of the following diagrams: 
We have two bulk diagrams presented on Figure 1: 
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where the bulk interaction point, denoted by an empty circle, represents z = ((, r) and we have 
to integrate over the whole left/right space-time. Thus the contribution of the first diagram 
is 



m 2 b 2 ,00 ,0 



/oo ru 
dr / dzGZ(y,z)GZ(z,z)Gt(z,y') 
-00 J —00 



Clearly G_(z,z) is divergent and we have to regularize it by introducing a cutoff A: 

f d 2 q i f A 1 dk . . 

G_{z,z)= / 777^0^ 2 , . = / ; 7T = A ( m cl) 

J (2vr) 2 g2 - ml + ie y ^ fc2 ~ m 2 2vr 

Its contribution can be absorbed into the renormalization of the mass parameter m\ 
m cl — m clT^( mcl ) wm ch results in extra counter term diagrams presented on Figure 2: 




Figure 2: Bulk counter-term diagrams with prefactor — m^A^i) 

The contributions from the defect terms can be grouped in two sets of diagrams. The first 
contains the same divergent loop integral and consists of those on Figure 3: 
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o 




o 



Figure 3: Divergent even defect loop diagrams with prefactor mc 8 ' b cosh \x 
where the interaction vertex is even together with the odd diagrams presented on Figure 4: 
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o o, 
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Figure 4: Divergent odd defect loop diagrams with prefactor m<: 8 1& sinh \i 



We have to integrate in time r over the real axis and the left / right part of the defect represents 
the contraction with the operators and 5> + , respectively. They all contain the divergent 
and regularized A(m c i) loop integral which can be absorbed into the renormalization of the 
defect parameter m c \ — ► m c \ — m c i^-A(m c i), which is consistent with the bulk renormalization. 
The resulting counter-terms produce the diagrams on Figure 5: 




Figure 5: Defect counter terms with prefactors — m c \\ A(m c \) cosh [i and — m c \\A(m c \) sinh /i, 
respectively 

The fact that all these singularities can be absorbed into the renormalization of m c \ is a 
nontrivial statement, since we have eight divergent diagrams having different propagators on 
the outer legs those we canceled just by renormalizing one single parameter in the original 
Lagrangian. The form of the renormalized Lagrangian is the same as the original one thus 
the integrable/topological nature of the defect is not spoiled by quantum effects, there is no 
anomaly. Observe also that the bulk mass term vr? cX and the boundary term m c \ renormalizes 
the same way so the bulk is the square of the other. 

The last group of the diagrams is the one which really contributes to the transmission 
factor. They are presented on Figure 6: 





Figure 6: Defect diagrams contributing to the transmission factor. They have prefactors 



■ sinh /j, and 



■ cosh respectively 
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The contribution of the first diagram is 



m c \b 2 



sinh fj, 



L 



oo 



dTGZ(y,z)G±(z,z)G±(z,y') 



4 



— oo 



Each of the terms on Figure 6 contains the finite contribution of the propagator 



Gl(z,z) 



I 



d 2 q i 



-T±(u,k) 



JL 

2vr 



(2-7r) 2 q 2 — m^j + ie 



This term together with the prefactor can be interpreted as the finite renormalization of the 
parameter — ► fj, + Sfi, where 6/j, = —^f- Summing up all the contributions and taking into 
account the different transmission factor dependent contributions on the outer legs we obtain 
the correction of order b 2 to the transmission factor as 



which is in complete agreement with (|23h . Thus we confirmed the renormalization of the 
parameter fj,, but we have not checked the renormalization of the parameter M which has not 
shown up at this order. We suspect, however, that in this perturbative scheme the boundary 
parameter m c \ renormalizes as the square-root of the term and only \i renormalizes as 
Bfi. It would be interesting to perform a two-loop perturbative calculation to decide about 
the renormalization of M c \ in the perturbative scheme. 

We note that we have also performed a perturbative calculation of order b 2 of the reflection 
factor, which can be extracted from G+(y,y'), and confirmed the absence of reflection at this 
order. In [39] the form of the sine-Gordon Lagrangian was fixed at order 6 6 by demanding 
the absence of particle creation. Following a similar line it would be tempting to see how the 
absence of reflection restricts the form of the defect potential in perturbation theory. 

5 Defect thermodynamic Bethe ansatz 

In this section we would like to check the consistency between the transmission factors and 
defect energies. In doing so we derive a DTBA to describe the ground state energy of a 
purely transmitting diagonal integrable defect on the circle of perimeter L. In boundary and 
defect systems there are two inequivalent ways to derive TBA equations. We can either focus 
on the groundstate energy or on the so-called ^-factors which is related to the finite volume 
normalization of defect/boundary states. Both BTBA equations have been analyzed in [To] 
although the (^-function type required further refinement |40| . Some details of the (^-function 
type DTBA can be found in [UJ and references therein. Here we consider the groundstate 
DTBA in the periodic setting as opposed to the strip geometry analyzed in |21j . 

5.1 Derivation of DTBA 

In order to derive the DTBA equation for the ground state energy we compactify the time- 
like direction with period R and calculate the partition function in two inequivalent ways by 
changing the role of the space and time coordinates. 

In the original description the defect is located in space as drawn in Figure 7: 



T-(9,B0(jj)) = T-(9, t i) + 



lib 2 1 

87r 1 - i sinh(6» + (i) 



+ 0(b A ) 
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i defect line 











Figure 7: Defect is located in space 



By taking the R —* oo limit the groundstate energy can be extracted from the partition 
function as 

lim Z(L, R) = lim Tr ( e- H ^ R ) = e~ E ^ L)R + ... 
In the alternative description when the role of time and space is exchanged as shown on Figure 



/L \ 



defect operator 



Figure 8: Defect is located in time acting as a defect operator 



the defect becomes an operator of the form [21] 

D = exp | J°° - 9)a + (9)a(9) 

which acts on the Hilbert space of the periodic model, TL. The partition function can be 
calculated as 

Z(L, R) = Tr (e-^D) = £ (24) 

\n)en 

The Hilbert space consists of multi-particle states 



(n\n) 



h > & 2 > ■ ■ ■ > Vn 
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on which the defect operator collects nontrivial diagonal matrix elements from 



?7T 7 7T 7 7T 

d|0i,0 2j . . . ,e n ) = t + {- - 0i)t+(- - e 2 ) . . . T + ( y - e n )\e 1 ,e 2 , ...,e n ) + ... 

while the energy operator acts as 

H\0i,0 2 , . . . ,9 n ) = (mcosh(^i) + mcosh(6> 2 ) + . . . + mcosh(6> n )) \9\, 9 2 , ■ ■ ■ , 9 n ) 

We can introduce another energy operator via H = H — \ log D such that the partition 
function can be written as 



Z(L,R) = Tr (e~ H ^ L D) = Tr (e'^ 1 ) = £ 



e 



-E„(R)L 



\n)eH 



This partition function can be calculated in the R — > oo limit by standard saddle point 
approximation taking into account the scattering of the particles. The calculation follows the 
usual route of TBA calculations, but now the kinetic term is shifted mcosh# — ► mcosh^ — 
-^logT + (^ — 9). As a consequence we obtain the following DTBA equations for the pseudo 
energy 

e(9) =mL coshO -log T+(— -6) - f — (f>(9 - 9') log(l + e'^) (25) 

2 J-OO ^ 

where <p(9) = —i4slogS(9). Once e is known the ground state energy can be expressed as 

/OO JQ 
— cosh# log(l + e~ m ) 
-oo 27r 

Here we do not have to shift the cosh 9 term since it comes from the derivative of the momentum 
(sinh0), which appears in the quantization condition. For the result in this generality see e.g. 

Alternatively, we can redefine the pseudo energy as e{9) = e(ff) — logT + (^ — 9) to obtain 



— 4>{9-9')\og[l + T + C—-9')e 

-oo 27 *" V 2 

from which the ground state energy turns out to be 

E (L) = -m ^ cosh^ log (l + T + (y - e)e~ £ ^ 

The ground state energy is real which can be easily seen form (|6|) since T£(^ — 9) = T + (^+9). 
As a simple consistency check we can see that the DTBA equation for the trivial defect, 
T + = 1, reduces to the periodic TBA equation [43]. We analyze the large and small volume 
limits separately in the next two subsections. 

5.2 Liischer type correction in defect systems 

If the volume L is large then e{9) = mLcosh# is large and we can expand the logarithm to 
obtain 

/oo f in 
— COsh^ T+( g^-mLcoshe + ( e -2mL ) 
-oo 2?r 2 
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This result can be calculated directly from (|24|) by taking the large L limit there. In that 
case, however, only the one-particle transmission term contributes which is universal for any 
quantum field theory. The result obtained is the analogue of the boundary Luscher type 
correction to the ground state energy [44] and is valid in any theory even in non-integrable 
ones. 



5.3 Defect energy- 
Here we analyze the small volume behavior of the ground state energy. Its normalization 
depends on the scheme in which the quantum field theory is defined. If we would like to 
compare the DTBA normalization to that of a perturbed defect conformal field theory, in 
which the perturbing operators have dimensions h on the defect and (h, h) in the bulk, then 
we have: 

2_ 00 

Eo(L) = -e de{ -e hulk L + T Y,c n l n{1 - h) ; I = mL 

n=0 

Only the perturbative terms, c n , are present in a perturbed rational defect CFT. (In non- 
rational CFT-s, like the UV limit of the boundary sinh-Gordon theory, we expect terms with 
logarithmic behaviour, see [Mj for the details). By calculating the small volume limit of 
Eq(L) from DTBA ed e f and et, u ik can be extracted exactly. The computation is analogous to 
the boundary one [HUGS], so we sketch only here. In the L — > limit the solution for e in ([25]) 
develops two kink regions around 9 = ± log | and a breather region around the origin. The 
behaviour of the solutions are determined by the 9 — > ±00 asymptotics of the integral kernel 
and defect source term: 

4>(9) = Ce~^ + 0(e~ 2 ^) ; log(T+(y - 6)) = A±e* e + 0{e^ 2B ) as 9 — ±00 

The two kink functions are responsible for the terms giving the central charge and the bulk 
energy constant, while the central/breather part gives the defect energy in the following form: 

m? m(A + + A-) 
etmik = '■> e def = 2c 

We note that the kink type behaviour does not exists for the whole parameter range of C 
and A±. The results is understood that we analytically continued it from a range where the 
calculation is reliable. 

Let us concretes the result for the two cases in question. In the sinh-Gordon model 

C = 4sinh7r.B ; A+ = -2e Tm 



so using (f26l) we recover {9} and JT8 
In the Lee- Yang case we have 



C=-2y/3 ; A± = tH^™^ + e ±i7r V) 



Plugging these expressions back to (1261) the results confirms the bulk energy density and the 
defect energy. 

We emphasize that the agreement obtained in the two cases confirm the solutions on one 
side and the DTBA equation on the other. 
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The other perturbative coefficients c n can be calculated from the DTBA only numerically. 
In the Lee- Yang case, however, one can gain further analytical information. One has to define 
Y{9) = e~ e ^ and to show from (j2~Tl) that it satisfies the Lee- Yang ^-system relation 



Y (e--)Y{e + -) = i + Y{e) 

from which the Y(ff) = Y(9 + ^p) periodicity follows. Similarly to the boundary case this 
gives the exponent of the perturbative expansion to be | showing that the dimension of the 
perturbing operator is h = — |. 

6 Defect bound-states and bootstrap closure 

In this section we analyze the analytic structure of the transmission factors for the whole 
range of their parameters. Since the sinh-Gordon transmission factors (fTBT) never have poles 
in the physical strip we focus on the Lee- Yang model only. Recall that in our convention the 
physical strip of the transmission factors T^{9) are Qm9 £ [0, |-]. 

6.1 Pole analysis on the ground-state defect 

We analyze the pole structure of both 

T_ (9) = [b + 1] [b - 1] and T+{9) = [5-b] [-5 - b] 

as the function of the parameter 6, simultaneously. We note that by folding the theory to 
a boundary one (with two particles) we could analyze its bootstrap in the usual boundary 
formulation. Here, however, we present the results in the defect language since the diagrams 
are more clear-cut. (The whole procedure, going from the reflection factor to the defect 
transmission factor, can be interpreted as taking a sort of square root of the boundary theory. 
By closing the defect bootstrap we would like to show that such a theory is indeed sensible. 
Since on the I boundary there is no boundstate any pole of the reflection factor appears either 
in T_ or in T + so the bootstrap will be very similar to the boundary one |6j). 

In determining the fundamental range of the parameter b we can see that b — > b + 12 is 
a symmetry. Moreover b «-> 6 — b exchanges T_ <-> T + so we can restrict ourselves to the 
range b £ [—3,3]. We will see by analyzing the defect excited states that the fundamental 
range is even smaller, only b 6 [—3, 2], as in the boundary case, since at b = 2 the role of the 
ground-state and the first excited state is exchanged. 

The poles and zeros of the transmission factor T_ (9) are 

poles of T_ are at 9 = —i—(b ± 5) ; zeros of T_ are at 9 = —i—(b ± 1) 

6 6 

The analogous expressions for T + are 

poles of T + are at 9 = i— (b ± 1) ; zeros of T + are at 9 = i—(b ± 5) 

6 6 

They can be drawn as the function of the parameter b as shown on Figure 9: 
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Figure 9: Poles and zeros of T + and T_ as function of b . Solid lines correspond to poles, 
while dashed ones to zeroes. The dotted lines show the fundamental range. 



For b G [—1,2] there is a pole in the transmission factor T + at 6 = iu = i^(b + 1), for 
which we associate a defect boundstate and denote it by |1+). Its energy is tocos ^(6 + 1) 
and the corresponding excited transmission factor can be calculated from the defect bootstrap 
equation shown on Figure 10: 




Figure 10: Defect bootstrap equations 



t! 1+> (#) = T_(6)S(9 + iu) 

From the defect crossing symmetry (|6j) we can calculate rJ. 1+> (0) as 

r{ 1+> (#) = t! 1+> (ot -6)= T_(ivr - 6)S{m - 6 + iu) = T + (9)S(6 - iu) 

which is consistent with the other bootstrap equation where the second particle arrives from 
the right. The resulting transmission factors are 

ri 1+> (0) = [6+1] [6 + 3] ; T{ 1+ V) = [5 -6] [3 -6] 

They are related to the groundstate ones as T± + \b —y 4 — b,8) = T T (b,8). This symmetry 
together with the defect energies indicate that when b exceeds 2 the role of the ground-state 
and the excited state |1+) are exchanged. This confirms that the fundamental range is indeed 
b €[-3,2]. 

In the range b S [1,2] the transmission factor T + {9) has another pole at 6 = i^(b — 1) for 
which we associate the defect boundstate |2+). It has energy mcos ^ (b— 1) and transmission 
factor 

T^ +) (d)=T ± (9)S(eTi^(b-l)) 

b 

Now we turn to the pole analysis of excited defect states. 
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6.2 Pole analysis on the exited defect state |1+) 

The poles and zeros of the transmission factors on the state |1+) are indicated on Figure 11. 



K 



2 








b 



-3 



-2 



-1 







2 



3 



Figure 11: Poles and zeros on the |1+) defect excited state. Dotted line shows the range where 
the excited state |1+) exists 



The state exist in the b £ [—1, 2] domain so we have to explain the poles in this range only. 

The pole of labeled by 1 on Figure 11 is at the same location as the one which 

creates the excited state itself, namely at 9 = if (6 + 1) in the full range b £ [—1, 2]. It can be 
explained by the first of the defect Coleman-Thun diagrams on Figure 12 



Figure 12: Defect Coleman-Thun diagrams for the |1+) state. The angles are measured in 
units of l -f 



The pole of T+ } labeled by 2 on Fi gure 11 is at 9 = if (6 + 3) and can be explained 
by the second diagram on Figure 12. Observe that by applying the Cutkosky rules [2] we 
would obtain a pole of second order but the transmission factor T_ has a first order zero at 
= if (1 — b) which, in this way, reduces the order of the pole to one. 

The pole of t!_ 1+ ^ labeled by 3 on Figure 11 is at 9 = if (3 — b). In the range b E [0, 1] it 
can be explained by the third diagram on Figure 12. Since the transmission factor T_ on the 
ground state has a zero at 9 = if (1 — b) the order of the diagram is reduced to one again. In 
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order for the diagram to exist the particle has to travel towards the defect, that is 1 — b > 0. 
This explains the pole in the range b 6 [0, 1]. In the range b 6 [1, 2] the particle creates a defect 
boundstate which is nothing but |2+). This can be seen both from the energy of the excited 
state m cos f (b + 1) + m cos f (3 — b) = m cos f (b — 1) and from the transmission factor. If the 
left particle creates a defect boundstate at rapidity 9 = iu then the excited states transmission 
factors are T±{9) = T±(9)S(9 ± iu). Now we can see from the bulk bootstrap equation that 

Ti 1+) (9)S(9±£(3-b)) = T±(9)S(9^l(b+l))S(9±£(3-b)) = T ± {9)S(9^(b-l)) = T^ +) {9) 

that is the transmission factors also supports the identification. 

6.3 The pole analysis on the excited defect state |2+) 
The defect boundstate labeled by |2+) has transmission factor 

Ti 2+> (#) = [b- l][6 + l] 2 [6 + 3] ; T{ 2+> (0) = [3 - b][5 - b} 2 [7 - b] 
The singularity structure can be summarized as follows. 
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Figure 13: Singularity structure of the transmission factors on the |2+) state. Bold straight 
lines represent poles of second order. The relevant interval where the boundstate |2+) exists 
is indicated by dotted lines. 



The pole labeled by 5 on Figure 13 is in T + (9) at 9 = i^(b — 1) and can be explained in 
the full range b £ [1, 2] by the first diagram on Figure 12, if we replace |1+) by |2+). 

124- *) 

The pole labeled by 6 on Figure 13 is in T_ (9) at 9 = if (3 — b) and can be explained 
by a diagram similar to the third one of Figure 12 in which the |1+) state is replaced by |2+) 
and the vacuum |0) is replaced by |1+). 
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Figure 14: Defect Coleman-Thun diagrams for the excited state |2+) 



The pole labeled by 7 on Figure 13 is a second order one in T + (9) at 9 = if (6 + 1) and 
can be explained by the two diagrams on Figure 14. Clearly the transmission factor T-{6) 
does not have zeros neither at 6 = if (3— b) nor at 9 = if (b — 1) so the pole is of second order. 

By now we explained all the poles of all the transmission factors of the ground and excited 
defect states. We used either the creation of a new defect boundstate or presented the appro- 
priate defect Coleman-Thun diagram which was responsible for the singularity. By finishing 
this procedure the spectrum become complete and we managed to define a sensible defect the- 
ory. It would be nice to check these findings by the defect truncated conformal space approach 
(TCSA). 

7 Conclusions 

We have demonstrated how the fusion idea can be used to solve topological defects in the 
sinh-Gordon and Lee- Yang models. In the sinh-Gordon case we determined the transmission 
factors and the defect energy as a function of a bootstrap parameter whose relation to the 
Lagrangian was also given. We checked these results in perturbation theory and against the 
newly derived DTBA. 

In the Lee- Yang case we determined the transmission factors together with the defect 
energy and checked them in DTBA. For certain range of the parameter the transmission 
factor admits poles in the physical strip. We closed the defect bootstrap programme: we 
explained all poles either by associating new defect boundstates or by giving the appropriate 
defect Coleman-Thun mechanism both for the groundstate and for excited defect states. 

The relation obtained between the transmission parameter and that of the Lagrangian in 
the sinh-Gordon theory can be analytically continued to describe the analogues relation in 
the sine-Gordon theory. This result also passes the test of first order perturbation theory 
and together with the transmission factors obtained in [22l E3J gives the complete solution of 
defect sine-Gordon model. We have checked this solution by performing the fusing procedure 
on the solitonic transmission factors. This is analogous to the dressing procedure in the XXZ 
spin chain developed in [46J. 

The perturbation theory developed here can also be used in higher rank affine Toda theories 
to connect the parameters of the transmission factor of the bootstrap solution [24] to the 
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parameters of their Lagrangians [18J. 

The derivation of the DTBA generalizes to any diagonal scattering theory. A large and 
small volume analysis analogous to the one presented in the paper will provide the leading 
finite size correction to the groundstate energy and give the bulk/defect energies, respectively. 

In the present paper we were concerned with the bootstrap (IR) description of our models. 
There is a need, however, to understand their UV behavior which probably can be described 
by perturbed defect CFTs. To connect these alternative descriptions we can use methods 
starting either from the IR side, like DTBA, or starting from the UV side, like defect TCSA. 
There are works in progress in both directions. 
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